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Abstract: In this paper we develop a method introduced by one of us to study 
metastable states in spin glasses. We consider a 'potential function' defined as 
the free energy of a system at a given temperature T constrained to have a fixed 
overlap with a reference configuration of equilibrium at temperature T' . We ap- 
ply the method to the spherical p-spin glass and to some generalization of this 
model in the range of temperatures between the dynamic and the static transition. 
The analysis suggests a correspondence among local minima of the potential and 
metastable states. This correspondence is confirmed studying the relaxation dynam- 
ics at temperature T of a system starting from an initial configuration equilibrated 
at a different temperature T' . 
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1 Introduction 



The off-equilibrium dynamics of glassy system is a fascinating subject. Experimental 
[p]] and numerical 0] evidence show that out of equilibrium phenomena persist in 
glasses for the largest reachable observation times. 

Recently a partial comprehension of these phenomena in mean field has been 
achieved || [|, || . The phenomenon of aging is described in the mean field theory as 
an asymptotic stationary state, where time translation invariance and the fluctuation 
dissipation relation do not hold. This regime turns out to be closely related to the 
nature of the static glassy transition. The models studied up to now can be divided in 
two classes according to their pattern of replica symmetry breaking (RSB) If the 
RSB is "continuous", i.e. if the Parisi order parameter function q(x) is continuous, 
then the static and the dynamic transitions occur at the same temperature. The 
asymptotic state is such that expected values of quantities which depend only on 
the configuration of the system at a single time (e.g. the energy or the distribution 
of the magnetizations) tend to their Boltzmann-Gibbs values. If instead the Parisi 
function is discontinuous, e.g. a single step function, the dynamic transition occurs 
at a temperature higer then the static one, and the quantities mentioned above tend 
to limits different to their canonical averages. However, a careful analysis in the 
spherical p-spin model || have shown how this limiting values can be related to 
the values of the same observables in a particular class of metastable solutions of 
the TAP equations, with free energy higher than the ones dominating the partition 
function. It then raises spontaneously the question whether it is in general possible 
to relate the properties of the asymptotic dynamic state to some static properties 
of the systems. 

This paper constitute a step forward in this direction. Our basic tool of in- 
vestigation will be a "potential function" introduced in j7|, defined as the minimal 
work needed to keep a system at temperature T at a fixed overlap from a typical 
equilibrium configuration of the same system at a different temperature T' . The 
models we choose to study are a family of spherical spin glass models which present 
a one step RSB. The most studied representant of this class is the spherical p-spin 
model. This model has the remarkable (and atypical) property that the order in free 
energy of the solutions of the TAP equations do not vary with the temperature. It is 
interesting to consider more general models where the order of the solution depends 
of T, and in the spin glass phase the low-lying states at different temperature are 
not correlated. 
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The basic theses of this paper is that the minima of the potential are related 
to (meta) stable states. To confirm this point we study the dynamics at tempera- 
ture T starting from a configuration which is perfectly thermalized at a different 
temperature T' at the initial time. 

Work on subjects related to the ones treated in this paper has been very recently 
achieved by R.Monasson ||. We thank him to make us aware of some of his results 
prior to publication. Applications of 'potential' model to Ising spin has also been 
done by 0. 

The organization of the paper is the following: In section 2 we present a short 
summary of the static and dynamic properties of the model. This chapter is part 
review and part new elaboration. Then in section 3 we introduce and discuss some 
basic properties of the potential function. Section 4 is devoted to the study of the 
potential for the spherical model. In section 5 we present the dynamic theory of the 
evolution of a system at temperature T starting from an equilibrium configuration 
at temperature T'. We finally sketch our conclusions. 



2 The Model 



The spherical p-spin model [10| is defined by the Hamiltonian 



H p [s] = Jh,...,i P s h-- s i P (!) 

n<«2<-<ip 

where the "spins" Si, i = 1,...,N are real variables subjected to the constraint 
(1/-/V) Yh s 1 — 1; an d the couplings Ji lt ... t i p are independent Gaussian variables with 
variance: J}, 



It can be observed that such model is (dynamically and statically) equivalent to 
a model in which the functional form of the Hamiltonian is not specified explicitly 
and one just take it as a random Gaussian function of the spin configuration with 
correlation function given by 

I 



H[s]H[s'] = N~q p StS , (2) 
g s , s ' = s ■ s'/N = (1/-/V) Yi s i s i is the overlap among the two configurations s and s'. 



One can then easily generalize the spherical model to random Hamiltonian Jy], |12 
with arbitrary correlation functions: 



H[s]H[s'] = Nf(q s>s ,). (3) 
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Polynomial functions f(q) with positive coefficients correspond to Hamiltonians rep- 
resentable as a sums of independent monomials of the kind ([!]). Depending on the 
function / the model can undergo either 'continuous replica symmetry breaking' 
or 'discontinuous' one. If the function g(q) = (f"(q))~ 3 ^ 2 f"'(q) is monotonically 
increasing with q in the interval [0,1], we are in the first case and we find a low tem- 
perature phase characterized by a continuous Parisi function If otherwise the 
function g monotonically decreasing one finds q(x) to be a single step function in the 
spin glass phase. Throughout this paper, if not stated otherwise, we will consider 
this second case; we will also always consider functions / such that /'(0) = 0. Here 
the asymptotic off-equilibrium dynamics and the statics lead to different results as 
far as the expect ed values mentioned in the introduction are concerned. Notably 
it is found for the energyat temperatures less than a dinamical critical temperature 
T D : lim^oo E Dyn (t) > E Gibbs , where E Dyn {t) is the energy computed in the infinite 
volume limit starting from a random initial configuration . 

For further reference we give here the result of standard static and dynamic 
analysis for the generalized spherical model. 



2.1 Statics 

From the standard study of the equilibrium measure of the model by means of the 
replica method it is find a free energy functional given (in standard notations) by 
the one step RSB form 

F = -flogZ 

= ~\ {/?[/ (1) - (1 - *)/(?)] + riog(l - Q) - t\ log [ 1 _\~l x)q ) } (4) 

q and x are variational parameters with respect to which F has to be maximized. 
In the high temperature phase q is equal to zero, and it is discontinuous at the 
transition point T = Tg, where x = 1 and q verifies: 

0sf(Q) + Ml - ?) + q = 0* s f\q) --*- = 0. (5) 

1 - q 



1 It has been observed in [[L2| that the model with f(q) = (1/2) (q 2 + q 4 ) has the same critical 
behaviour the SK spin-glass. 
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2.2 Off-equilibrium dynamics 

One studies here the Langevin dynamics associated to the model with random initial 
condition. The system is analyzed supposing that the thermodynamic limit is taken 
before the infinite time limit, so as to prevent full equilibration in presence of a phase 
transition. 

The relevant objects of investigation (order parameters) turn out to be the cor- 
relation function C(t,t') = l/NJ2i s i(t) s i(t')i an d its associated response function 
G(t,t') = l/NJ2i th-ft 1 ) ( nere an d in what follows t > t'). One finds then the set of 
coupled equations: 

-n(t)G(t, t') + f ds f"(C(t, s))G(t, s)(G(t, t') - G{s, t')) , (6) 
Jo 

-nC(t,t') + C ds f'(C(t,s)) G(t',s) 
Jo 

fds f"(C(t,s))G(t,s) (C(t,t')-C(s,t')) , 



fds f'(C(t,s)) G(t,s) 



r ds f'\C(t, s))G(t, s) (C(t, t) - C(s, t)) + T. (7) 
o 

Following the analysis of || for the p-spin model, one finds that for large times, 
two regimes are important: a first one in which one takes the limit t,t' — ► oo fixing 
to a finite value the difference r = t — t' and a second one where the limit is taken 
fixing the ratio A = h(t')/h(t). h is a function that at present the theory is not able 
to specify. 

• In the first regime the correlation function decays with r from the value one 
to a finite value q (dynamic Edward- Anderson parameter) and the fluctuation 
dissipation relation TG as (r) = —dC as /dr is satisfied. 

• In the second regime one finds a non time homogeneous form C(t,t') = 
C(h(t')/h(t)), with C(A) monotonically increasing with A and C(l) = q and 
C(0) = 0. The response function is equal in this regime to TG(t, t') = x 9C ^f\ 
The constant x is found to be a number between zero and one in the low 
temperature phase. The parameters q and x are solution of the equations: 

P 2 f"(q)(l-q) 2 = l. (8) 



dG{t, f) 
dC(t,t') 



fi(t) 



+ 



+ 



x 



qf"(q) 
/'(?) 
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The dynamical transition is marked by the condition x = 1 at T = Tp, and 
one can check that Tp > T$. The first of eq. (§) is equivalent to the variational 
equation dF/dq = in statics, while the second one coincides with the condition 
that 'replicon eigenvalue' of the fluctuation matrix in replica space is equal to zero. 
We will refer to that equation as marginality condition. 



2.3 TAP approach 

For completeness we mention a third approach that is useful to investigate the sys- 
tem. This is the one of the TAP equations, in which one writes mean field equations 
for the magnetizations for fixed disorder. Using the diagrammatic approach of | ill, 



based on the skeleton expansion introduced in fll4| , one easily finds the TAP free 
energy function: 

Ftap [ m > q] = ±H[m] - ~log(l - q) - - f{q) - (1 - q)f(q)} (9) 

where the variables m 8 represent the average magnetizations, and the self-overlap 
q is given by l/NJ2i m i- The physical states are solutions of variational equations 
both with respect to the m's and q. 



In the p-spin model it was found useful |L5] to rescale the variables to 

m -> ^fqsi -^I>i=i- (10) 

i 

It then follows from (|1|) that 

H[m] = q p/2 H[s] (11) 

The points of extreme of Ftap with respect to the S{, as well as their order in free 
energy do not depend on temperature |1J. Let us stress here that this is a highly 



non generic property which depends critically on the homogenity of the Hamiltonian 
(0). In general we can expect the order of the solutions to depend on temperature. 
The situation in which this happens for the lowest states is usually called chaotic 
in the literature fl6f , and has been recently fully demonstrated in the SK model 



A complete analysis of the TAP equations equations is at present missing. In 
section 4 we will find evidence for chaos with respect to temperature in spherical 
models with the method we propose in the next section. 
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3 The Method 



We now introduce the "potential" , our basic tool of analysis in this paper. 

We consider an arbitrary configuration of the spins s, drawn as a typical re- 
alization of the canonical probability distribution at temperature T', Pcan[s] = 
exp(—/3'H[s])/Z(T'). WE can then compute the cost in free energy at a tempera- 
ture T (in general different from T') to keep the system at a fixed overlap ft = q s a 
with s; namely 

V = -^logZ[s;p}-F[T}; (12) 
Z[a;p] = J daexp(-(3H[a])5(p-q Sta ) (13) 

where F[T] is the free energy without constraint. As Z in (|13"D is a sum of positive 
terms, V in (|T2]) is a positive quantity. It is reasonable to suppose that V in addition 
of being self- averaging with respect to the quenched disorder in the Hamltonian, it 
is also self-averaging with respect to the probability distribution of the reference 
configuration s. In this problem the spins Si are quenched variables on the same 
foot as the random Hamiltonian itself. In this respect the present method allows to 



extend and improve the analysis of ref. ||18|| , where there were considered different 
'real replicas' coupled in a symmetric way. 

We need then, in order to compute V, to perform the average over the Hamilto- 
nian and the reference configuration s. 



NV = ^ J dsexp(-P'H[s])[-T\ogZ[s;p] - F[T]]. (14) 

This average can be done with the aid of replica method. We found two strategies 
which consistently lead to the same results. The first one is based on the formula 



NV = -T\imhm — ^exp(-P'H[s])Z[s-,p]R-A . 



(15) 



One then evaluates the average for integer R and n and obtains V by analytic 
continuation. 

The starting point of the second strategy is the formula 



NV = -T lim lim V exp (-(3'H\s}) Zip'}™- 1 \ Z ^^ m 1 



(16) 



and again one performs an analytic continuation from integer n and m. The use of 
one procedure instead of the other is mainly a matter of taste. We will sketch the 
first stages of the formal manipulations for both the procedures, and we will treat 
our model with the second one. 

Let us start with the first procedure. The replicated partition function is: 



z (n,R) = J dgajr 



TL Yl R 

a=l a=l r=2 



a=l r=2 



Npj. 

(17) 

One can then perform the average over the distribution of the Hamiltonian, and 
introduce the order parameters: 



Q x \ ^ r,a s,b 
(o,r),(M) = Jjl^ S i S i ■ 



(18) 



The logarithm of Z™ ,R divided by N is then found to be equal to the saddle point 
over the Q's of 



^EEW(QM,(M)) + ^logQ 

Z r,s a ,b Z 



(19) 



where (3 r = (3' for r = 1, f3 r = (3 for r = 2, ...,n. As usual one needs a scheme to 
perform the analytic continuation in the number of replicas. 
With the second procedure, one has 
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ds a J da a exp 



a=l a=l 



( 2 °) 

a=l V i / 



After the average over the distribution of the Hamiltonian is performed, one 
introduce the order parameter matrices: 



1 



Qab = j^J2 s-s\ 

i 

R *I3 = a i a i 



1 



N ^ 



(21) 
(22) 
(23) 



with a,b — 1, n and a, (3 — 1, m. Combining the order parameters in the single 
[n + m) x (n + m) matrix 

'Q p\ 

P T R 



Q 



(24) 



one finds 



N l ° gZ2 "2 



+ iTrlogQ. 



l,n l,n l,n 

£/3' 2 f(Q a , b ) +J2/3 2 f(R a>p ) + 2j2{3/3'f(P a , a ) 

a,b a,/3 a,a 

(25) 

Observing that the constraint implies Pi )Q = p for any a = l,...,m, a sen- 
sible ansatz is to assume that the matrix P has elements P a , a = Pa indepen- 
dent of a for any a. A simple computation based on the formula Trlog(Q) ~ 
log(/ dxexp(— xQx)) reveals then the identity: 

Tr log[Q] = Tr log[Q] + Tr log[i? - A] (26) 

where A is a m x m matrix with all the elements equal to A a> p = J2 a b Pa{Q~ 1 )abPb- 
Note that with this ansatz for P, TrlogfQ] is of order n, while Trlogfi? — A] is 
of order m. So, neglecting terms of order m one finds that the equation specifying 
the matrix Q is just 

P' 2 f'(Q ab ) + (Q-% = (27) 

which is independent of P and R and is just the saddle point equation for a system 
at equilibrium at temperature T' . This is a good consistency check for the ansatz: 
the s system, which is at equilibrium, should not be affected by the a system. The 
variational equations for P and R are respectively written as 

BV 

-wZr = /W^)+£(Q-%A£((i*-^)-^ = o (28) 

b a,p 

BV 

"20 op- = i3pf , {R a p) + ({R-A)- 1 ) afl = (29) 

For a = 1 the l.h.s. of equation ( p8|) represents the derivative of V with respect to 
p and has not to be equated to zero. We note that for p = the solution to (p8|) is 
simply P a = for any a. One has then ^ = and the potential takes its minimal 
value V — 0. 

To solve the equations (^,^) for n, m — > 0, it is of course needed a continuation 
scheme for the various order parameter matrices. For Q the usual hierarchical 
ansatz has to be employed. A sensible choice is to take also the matrix R 
of the hierarchical form (with associated Parisi function r(x)), and the vector with 
components P a , as the first line of a hierarchical matrix P* b associated to the diagonal 
element p and to a function p(x). Although this general form can be of interest for 
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many probems, notably whenever RSB is necessary to find the free energy, we will 
see later that a replica symmetric ansatz will be enough for the problem addressed 
in this paper, at least for the class of models we consider. 

In the next section we will devote a lot of attention to the points of minimum of 
the potential V(p). The basic theses of this paper is that the minima of the potential 
correspond to metastable states. 

Before leaving this section we just comment on the fact that in the case of 
continuous RSB transition where the probability distribution P(q) is different from 
zero in a whole interval q min < q < q ma x, we checked that as it has to be expected, 
the potential is flat and equal to zero in the whole interval [q m in, qmax}- 



4 The potential in the intermediate regime. 

We study in this section the properties of the potential when the system s is in the 
high temperature phase, and the order parameter matrix is replica symmetric and 
given by Q a>0 = 5 a ^. According to the picture proposed in [19|, |T3[ for T D < T 
this replica symmetric phase does not describe an ergodic phase, but a situation in 
which there are an exponentially large number of states, with zero overlap among 
each other ||19|| . The triviality of the Parisi probability function P(q) = S(q), implied 
by the replica symmetric solution, is just due to entropic reasons. It is then natural 
to suppose for P and R a replica symmetric structure with 

Ra,b = S a , b + (1 - 6a, b )r. (30) 
The insertion in ( ^5l) gives for the potential 



v= ~h \ m '^ p) ~ /?2/(r) + log(1 _ r) + t^~} (31) 

where r is the only parameter with respect to which the potential has to be maxi- 



mized. The maximization equations ( p8| , |29D reduce to the single equation 



^V) = (32) 

It turns out that for any concrete form of the function / that we have analyzed that 
this equation admits a unique solution in the range of j3' where the replica symmetry 
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is not broken. As we stressed in the previous sections the points of extremum of V 
with respect to p are of particular importance in our analysis, the minima correspond 
to stable or metastable states in which the system can be trapped. Let us then write 
the equation that specifies the points of extremum: 

0f'(P) = t~ (33) 
1 — r 



We stress again that the point r = p = is always a solution of (p3| , p2|) , and is 
always an absolute minimum with V = 0. 

Let us start our analysis considering the case of equal temperatures T' = T. In 
this case the qualitative features of the potential are largely model independent. In 
figure 1 we show the picture of the potential for four different temperatures in the 
case f(q) = q 3 /2, similar plots are obtained for arbitrary functions. From top to 
bottom, they represent the potential at temperature higher then T D , equal to T D 
between To and T$, and right at T$. We can see from the figure that for T > To the 
potential is monotonically increasing, and the only extremum of the potential is the 
minimum at p = 0. At the temperature To where the dynamical transition happens, 
the potential develops for the first time a minimum with p = r. It is interesting to 
observe that the energy in this flex point | is equal to the asymptotic value of the 
energy in the off equilibrium dynamics. The same is true for the parameter r which 
turns out to be equal to the dynamical Edward-Anderson parameter. 

The condition for the potential of having a flex coincides with the marginality 
condition. Indeed the flex implies a zero eigenvalue in the longitudinal sector and at 
x = 1 the replicon and the longitudinal eigenvalues are degenerate (see for example 
the formulae in ref. |21J). This marginality condition is well known to give exact 



results for the transition temperatures in p-spin spherical models. It also give ac- 
curate results, compared with the Monte Carlo, in the case of the Ising model with 
random orthogonal matrix ]20| (ROM model). 

We have observed that in general more then one minimum can be present in 
the potential. In the p-spin model it happens that two minima develop at the same 
point. The rightmost one, that we will call primary is the one with p = r, while the 
other, secondary, has p < r. For temperatures smaller than To the minima have a 
finite depth, i.e. are separated by extensive barriers from the absolute minimum. 

The primary minimum is easily interpreted. There the system denoted by s is 
in the same pure state as the system a. In the region T$ < T < To the number of 
2 See the following for the definition. 



11 



0.04 



0.035 



0.03 



0.025 



0.02 



0.015 



0.01 



0.005 



-0.005 




Figure 1: The potential function in the p = 3-spin model for T = V and various 
temperatures. From top to bottom (3 = 1.59, f3 — f3 D — 1.63,/? = 1.66, f3 — Ps — 
1.70,. Observe that a minimum is first present at the dynamical temperature, and 
the value of the minimum of potential at the statical temperature is zero. 
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pure states is exponentially large in N: J\f = e N ^ T \ Consequently the probability 
of finding two system in the same state is exponentially small and proportional to 
e -A r s(T)_ rj^g f ree ener gy cos i t constrain two systems to be in the same state is 

then proportional to the logarithm of this probability, namely we have 

^primary THu(Tj. (34) 

Coherently at the statical transition temperature T = T$ one finds V pr i mary = 0. 
The quantity £ has been computed for the p-spin model in ref. [T3j as the number 



of solution of the TAP equation with given free energy and coincides with our 
calculation. The secondary minima, could aso be associated to metastable states, 
but at present we do not have an interpretation for them. This conclusion on the 
equivalence of the potential with the number of solution of the TAP equation hold 
also in the ROM J2TJ, and it has been argued by in ref. || on general grounds. 



The study of the potential for temperatures smaller than T5 would require to 
take into account RSB effects, which would complicate a bit the analysis. However 
it is physically clear that the shape of the potential in that region it is not different 
qualitatively from the one at T = Tg. It has a minimum where r = p are equal to 
the Edwards Anderson parameter and the value of potential is zero. 

In the case of different temperatures T ^ T' the properties of the potential 
depend on the presence or absence of chaos with respect to temperature changes. 
The primary minimum of the potential if it exists, reflects the properties of the states 
which are of equilibrium at temperature T' when they are followed at temperature T. 
In the p-spin model, where the order of the levels do not depend on the temperature, 
and each level has temperature independent complexity, the value of the potential 
in the primary minimum can be related to the properties of the solutions to the 
TAP equations when they are followed in temperature. Denoting by Ftap{T, E') 
the free energy of the TAP states that dominate at temperature T' when they are 
followed at temperature T, we find: 

Vprimary = -TE(T') + F TAP (T, E') - F{T) (35) 

where F(T) is the free energy at temperature T and Ftap{T, E') is the TAP free en- 
ergy at temperature T of the states which are of equilibrium at temperature T". We 
also note that q is equal to the EA parameter of the aforementioned TAP solution. Of 
particular interest is the case T' = Tp where the primary minimum is marginally sta- 
ble for any value of T, and has energy, defined as E prirnary = (d/dP)(3[V primar y—F(T)} 
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Figure 2: The potential for the p = 3-spin model and different temperature. The 
curves are drown for (3 = (3 S and, from top to bottom, (3' = 1.59, (3' = (3d = 1.63, f3' = 
1.66,(3' = (3 S = 1-70. 

and EA parameter equal to the these of the off-equilibrium asymptotic state. In fig- 
ure 2 we show the potential for p = 3 at fixed (3 = ((3 D + (3s) /2, and various values 
(3'. For a given T' the minimum exists in the range of temperatures for which the 
corresponding solutions to the TAP equations exist. 

In the case models with inhomogeneous Hamiltonian the situation for T ^ T' is 
different. We did not try a systematic study of this case and we concentrated on 
the case of a 3 + 4 spin model with f(q) = (1/2) (q 3 + eg 4 ). The first thing that we 
note in this case is that the horizontal flex present for T — T' — Td disappears as 
soon as, for T' = Td we take T ^ T/> For T' < Td it exists an interval of values of 
T around T' in which the minimum exists. In figure 3 we display the potential of 
the 3+4 model for fixed T' < T D and various values of T. 

One might wonder if the point of horizontal flex of the potential are associated 
to the dynamical state as it happens in the p-spin model. To do that we have looked 
at values of T and T 1 such that the marginality condition (3 2 (1 — r) 2 f"(r) = was 
verified in the minimum of the potential. In that points we found that the minimum 
has a finite depth and the energy does not coincide with that of the dynamical 
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states. 

It is tempting to interpret the minima of the potential as the properties of equi- 
librium TAP solutions at temperature V when followed at temperature T, we are 
however very cautious on this point due to the possibility of crossing of the so- 
lutions in free-energy. Let us just mention that this should be valid below the 
static transition (T" < T s ) where the lowest TAP states dominates the partition 
sum, and where there can not be crossing. Chaos with respect to temperature 
implies that the free energy of these solutions at temperature T will be typically 
higher then that of the equilibrium solution. One could expect that the minimum 
of the potential has a value higher then zero for any T ^ T' . We aspect then that 
V primar y(T,T') « (T - T') 2 \ with integer k. 



5 Dynamics 

We study now the relaxation dynamics at temperature T of a system starting its 
evolution at time zero from an equilibrium configuration at temperature T'. We 
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start from the usual Langevin equation for the model 

dsAt) , . , . OH , . 

— = -i I (t) Si (t)- — + Vi (t) 

(Vi(t)Vj(t')) = 2T5 ij 5(t-t') 



(36) 



As usual the time dependence of fi(t) is chosen to enforce the spherical constraint 
|s| 2 = N at any time. The dynamical generator of the correlation function [22], 
averaged over the initial condition is 

ZDyn 



1 



J V V VsVs^^ex V {-P'H[s(0)]}- 



exp 



/ dt' E isi(t) [si{t>) -\- — h fi(t)si(f) + TiSi(t')] 

JO ; OSi 



(37) 



The term 1/Z(T') obliges us to introduce replicas to perform the quenched average. 
One possible way is to start from the average of the logarithm of Z^ yn and replicate 
the system at each time |[23|| . We follow here a different route, using the relation 



1/Z = lim^o^™ -1 and replicating only the system at time zero. So, "replica 
number 1" will be present at all times, while replicas 2 to n will only be present 
at time zero. We will denote as sf(0) the spins of the a-th replica at time zero. 
Coherence would require the notation s\(t) for the value of the spins at time t; we 
will use instead Si(t). Using standard manipulations to average over the disorder we 
find 



J Dyn 



l,n 



E^/Wo)-^)/^)-! E*?(o) s 

a,b a,i 



VsVs exp 
exp f dt' E §i (f) [is, (*') + Ti Si (t')] 

exp [ t dt'dt"f , (s(t') ■ s(t")/N)Yis i (t')is i (t")+ 
l Jo i 

f"(s(t') ■ s(O/JV)s(0 • zs(t")/Nj2^(t') Sl (t" 

i 

exp U'E fdt'f{s{t>) ■ s o (0)/A0E^(*'K(0) 



(3? 



Introducing the correlation function C(t', t") = s(t')-s(t")/N, the response G(t' , t") 
s(t') ■ is(t")/N (t' > t") and the correlation with the initial condition C a {t') 
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s(t') ■ s a (0)/N, and observing that they are non-fluctuating quantities for iV — > oo, 
we recognize in (38) the dynamical generator functional of a system of N indepen- 
dent spins subject to the Langevin equation 

ds(t) r l n 

-Ai = -fi(t) Si (t) + / f'(C(t, t'))G(t, t') Si (t') +(3'J2 /'(C«(*)K(0) + Bi(t) 
dt Jo a=1 

(BityBjif)) = 2T5 il 5(t - + 5 h3 f'(C(t,t')) (39) 

satisfying the self consistency equations: 

C(t,t') = s{t)-s{t')/N 
G(t,t') =s(t) -is(t')/N t>t' 

C a (t) = s(t) ■ s a (0)/N. (40) 

It is easy to check that the "equilibrium part" in the action, involving the parameters 
Qab = s a(0) ■ s b (0)/N is not affected by the dynamics and, as it should, describes 
correctly the equilibrium at temperature T'. From the Langevin equation one derives 
the closed system: 



-fi(t)G(t, t') + f ds f"(C(t, s))G(t, s)(G(t, t') - G{s, t')) 
Jo 



dG{t, f) 
dt 

= -n(t)C(t, + f ds f(C(t, s)) G(t', s) 
at Jo 



+ fds f"(C(t,s))G(t,s) (C^O-C^O) +(3'Y,f(C a (t))C a (t / ) 

-n(t)C a (t) + f ds f"(C(t, s))G(t, s)C a (s) f'(C b (t))Q ajb 
Jo i 

fi(t) = [ ds f'(C(t,s)) G(t,s)+ [ ds f"{C{t,s))G{t,s) (C(t,t) - C(s,t)) 



dCg{t) 

dt 



+T +(3'J2f'(C a (t))C a (t), (41) 

a 

with the conditions: 

C 1 (t) = C(t,0) C a (0) = Q ha (42) 

These equations differ from the usual off-equilibrium equations for the presence of 
terms which couple with the configuration at time zero. The matrix Q a>b has the 
usual hierarchical form, and C a (0) is its first row. It is then natural to assume that 
C a (t) continues at any time to have the structure of the first row of a hierarchical 
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matrix. Once this ansatz is plugged in (4~I) we get a system of equations which 
admits for any time a unique solution. 

The basic question that can be answered with the aid of these equations is 
whether starting from equilibrium at temperature T' the system equilibrates in a 
metastable state, or it ends up in an aging state. In this paper we will only consider 
non-aging solutions to the equations (flip . We leave to future work the investigation 
of possible aging solutions. 

We look then for solution where no aging appears, the correlation functions 
C(t,f) tends for infinite times to the homogeneous function C as {t — £'), and the 
response function is related to C as via the fluctuation dissipation relation: G as (r) = 
—(3 9C q 3 t ^ ■ In the high temperature regime C a (0) = and the equations coincide 
with the off equilibrium ones. Interesting phenomena appear for temperature T' < 
To- Let us study the case T$ <T' < To- Here we have Q a ^ = 6 a ,b, a bit of reflection 
on the equation for C a (t) reveals, that the replica structure of C a (t) is at any time 
of the kind 

C a (t)=5 1>a C(t,0). (43) 



Using then the notations: 



p = lim C(t, 0) 

t— >oo 

r = lim C as (r) 



/Xoo = lim n(t) (44) 



one finds for C as {r) the equation 

= -^C as {T)+(3[f{l)C as {T)-r{T)r]-(3 [ nC as (rW))^P-+(3'pfip)- 

(45) 

While p, r, /ioo are solution of the following system 

-fXooP + PPPooP' f {$) = 

-/W + Wf'W - f'(r)} + (3(1 - r)f'(r) + f3'pf'(P) = 

+ T + (3[f'(l) - f\r)r] + P'pf'iP) = 0. (46) 

Eliminating jj,^ we find the remarkable result that the equations specifying r coincide 
with the equations ( |32| , |3"3"D of section 4 for the points of extremum of the potential. 



We have seen in section 4 that depending on the form of / and on the temper- 
atures T.T' solutions to (EBT) can exist or not. In the case in which such solutions 
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Figure 4: The convergence of C(t, 0) to p for the 3+4 model in the case T = 0.673, 
T' = 0.804. We plot C(t, 0) — p as a function of t in a log-log scale, for p we take 
the value in the primary minimum of the potential p = .646. C(t, 0) — p seems to 
approach zero as a power law. A similar plot is found for the energy as a function 
of time. 

do not exist, we easily conclude that for long times an asymptotic aging regime sets 
in the system. To understand if this regime is correlated with the initial condition 
would require a more complete ansatz then the one we are using here and we leave 
it f or future work. In the case solutions with p ^ exist, we need to show that 
the solution to (f|l]) converge to one of the solution of (|46|) and to know to which 
one. The problem is easily solved for equal temperatures T = T' . If we start at 
equilibrium we stay there. We explicitly verified that the solution to the equations 
(|4i"D is time translation invariant from the start. We belive that it a general fact that 
the dynamics tends toward the primary minimum of the potential. To verify that, 
we performed the numerical integration of the system fl4"I|) for f(q) = 1/2 (q 3 + q 4 ) 
with the method employed in ||. Figure 4 shows the approach of C(t, 0) to p = .646 
in the case T = 0.673, T' = 0.804. In the same way it is possible to see that the 
energy tends to the value predicted by the potential theory. 
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6 Conclusions 



In this paper we have investigated the free-energy landscape of a spherical spin-glass 
model, in the vicinity of equilibrium states. 

The model is chosen in a way to have a low-temperature glassy phase charac- 
terized by a single step Parisi function q(x). The study is performed introducing 
the 'potential function' defined as the minimal work needed to keep a system to a 
given distance from an equilibrium configuration. In this way we could identify the 
dynamical transition temperature as the point where first appear a minimum in the 
potential function at equal temperatures. 

In the p-spin model, where the Hamiltonian is an homogeneous function and there 
is not chaos with respect to temperature changes, the dynamical states for T < To 
can be associated to the points where the potential has a flex. This property is not 
true any longer whenever chaos is present in the model. In that case, we do not have 
a recipe to associate properties of the dynamic asymptotic state to some particular 
point in the phase space of our coupled system. 

The minima of the potential have been associated to metastable states. The 
life of the metastable states can be estimated arguing that the free-energy barrier is 
just given by the maximum to minimum difference in the potential 0, which is then 
proportional to the size of the system N. It is interesting to remind here the result 
of the analysis in [[?[] in finite dimension. In that case the order parameter is space 
dependent. The free energy barrier can be estimated with instantonic techniques 
(see e.g. ref. f24|), introduced in the the contest of spin glass models in p5| . It turns 
out that in three dimensions the relaxation time of the metastable states scales as 



Tmax ~ e( T - Ts ' 2 , to be compared with the Vogel-Fulcher law T max ~ e( T - Ts > , often 
used to fit the temperature dependence of the viscosity near the glassy transition in 
structural glasses. 

Let us conclude the paper noting that simple generalizations of the potential 
method would allow to gain further insight on the free-energy landscape of spin- 
glass models. For example one can consider the free energy cost of a situation where 
there is a first replica at equilibrium at a given temperature, a second one, at another 
temperature, is kept at a fixed distance from the first one, and a third replica, at a 
third temperature kept at fixed distances from the first two. 

Also, it would be very interesting to generalize our calculation to models with 
continuous RSB, to study the property of chaos with respect to temperature changes 
in that case. 
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